home *** CD-ROM | disk | FTP | other *** search
/ SGI Developer Toolbox 6.1 / SGI Developer Toolbox 6.1 - Disc 4.iso / src / exampleCode / MP / kernel_switch / main.f < prev    next >
Encoding:
Text File  |  1994-08-02  |  1.9 KB  |  96 lines
  1. c 
  2. c Copyright (C) 1993, Silicon Graphics, Inc.
  3. c All Rights Reserved.
  4. c
  5.       PROGRAM MAIN
  6.       IMPLICIT DOUBLE PRECISION(A-H, O-Z)
  7. C 
  8. C---------------------------------------------------------
  9. C 
  10. C       DRIVER PROGRAM FOR INNER PRODUCT ROUTINE
  11. C 
  12. C------------------------------------------------------------ 
  13. C 
  14. C 
  15.       DIMENSION A(800,800),B(800,800),C(800,800)
  16.       REAL*4 T0,ETIME,TARRY(2)
  17.     DOUBLE PRECISION ALPHA,BETA
  18.  
  19.     character*80 evalue
  20.     integer charcmp
  21.  
  22.     ALPHA=1.0D0
  23.     BETA = 0.0D0
  24. C 
  25. C 
  26.  300  READ(5,*,END=999) MA,NA,MB,NB,MC,NC,IDUMP
  27.       PRINT *,'MA,NA,MB,NB,MC,NC,IDUMP = ',MA,NA,MB,NB,MC,NC,IDUMP
  28. C 
  29. C 
  30. C       setup jump table
  31. C 
  32.     call readenv('SGI_ISA',evalue,7,length)
  33.     if (length .eq. 0) then
  34.        print *,'error - set SGI_ISA environment variable'
  35.        stop
  36.     end if
  37.     if (charcmp(evalue,'mips1') .eq. 0) then
  38.          call setmips1()
  39.     else if (charcmp(evalue,'mips2') .eq. 0) then
  40.         call setmips2()
  41.     end if
  42.  
  43. C
  44. C
  45. C       GENERATE A,B
  46. C
  47.  
  48.        DO 100 I = 1, MA 
  49.     DO 100 J=1,NA
  50.  100    A(I,J) = 1.0D0
  51.        DO 200 I = 1, MB 
  52.     DO 200 J=1,NB
  53.  200    B(I,J) = 1.0D0
  54. C       DO 300 I = 1, MC 
  55. C    DO 300 J=1,NC
  56. C 300    C(I,J) = 1.0D0
  57. C
  58. C       INITIALIZE NUMBER OF COUNTS
  59. C
  60.         ICOUNT = 2 * MA * NA * NC * IDUMP
  61. C
  62. C       FIND SYSTEM ID AND # OF HEADS
  63. C
  64. C
  65. C       CALL MATRIX PRODUCT ROUTINE
  66. C 
  67.       T0 = ETIME(TARRY)
  68.     DO III=1,IDUMP
  69.     CALL BLOC_MPY('N','N',MC,NC,MB,ALPHA,A,800,B,800,BETA,C,800,
  70.      +    1)    
  71.     ENDDO
  72.       T0 = ETIME(TARRY) - T0
  73.  
  74. C
  75.  
  76. C      OUTPUT RESULTS
  77. C
  78. C      WRITE(6,9) S1,S2,S3,S4,S5,S6,S7
  79. C 9    FORMAT(' ******** THE INNER PRODUCT ARE ',6G15.7/1X,G15.7)
  80. C      IF(IDUMP .GT. 0) CALL MATPRT(MC,NC,C,800,800)
  81.       CTIME = 1.0
  82.       WRITE(6,20) CTIME,T0 
  83.   20  FORMAT(//'* CPU TIME FOR OUTER PRODUCT =',G15.7
  84.      + /'  WALL TIME FOR OUTER PRODUCT =',G15.7)
  85. C
  86. C     CALCULATE MEGAFLOPS
  87. C
  88.       WFLOPS = DFLOAT(ICOUNT)/1.D6/T0
  89.       WRITE(6,31)  WFLOPS
  90.  31   FORMAT(' MEGAFLOPS =',G15.7)
  91. C
  92.       GO TO 300
  93. C
  94.   999 STOP
  95.       END 
  96.